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SintlARY 


The  development  and  the  numerical  application  are  presented  of  a 
Rayleigh-Rltz,  or  modal,  type  of  flutter  analysis  which  takes  into 
account  three-dimensional  structural  and  aerodynamic  behavior.  The 
flutter  mode  is  approximated  by  a  series  of  natural -vibration  modes, 
and  the  aerodynamic  forces  corresponding  to  these  modes  are  derived 
from  subsonic  lifting-surface  theory,  according  to  the  kernel-function 
approach,  for  a  finite  wing  oscillating  in  conpresslble  flow. 

The  application  is  made  to  a  delta  semispan  wing  with  a  leading- 
edge  sweep  angle  of  45°  which  fluttered  at  a  Mach  number  of  O.85.  Results 
of  flutter  calculations  show  that,  for  this  case,  when  the  first  three  or 
four  natural -vibration  modes  are  used  to  approximate  the  flutter  mode, 
converged  solutions  for  the  flutter  speed  are  obtained  that  are  about 
5  percent  less  than  the  experimental  value.  Theoretical  flutter-speed 
boundaries  were  located  for  a  range  of  densities  and  Mach  numbers 
including  those  of  the  experimental-flutter  condition.  Further  applica¬ 
tion  of  the  analysis  to  study  the  effects  of  variation  in  certain  struc¬ 
tural  properties  showed  that  the  converged  flutter  speeds  were  more  sen¬ 
sitive  to  variations  in  the  natiiral  frequencies  than  to  either  variations 
in  mass  or  to  the  inclusion  of  generalized-mass  coupling  terms  \diOse 
existence  is  due  to  the  use  of  e3q>erimental  natural  mode  shapes. 


INTHODUCnON 


Current  aircraft-design  trends,  such  as  the  use  of  thinner  struc¬ 
tures  and  external  stores  on  aircraft  capable  of  very  high  speeds,  have 
combined  to  diminish  flutter  safety  margins  and  have,  consequently, 
increased  the  need  for  greater  accuracy  In  flutter  prediction.  As  a 
result,  both  the  structural  and  the  aerodynamic  aspects  of  the  flutter 
problem  should  be  treated  by  use  of  more  realistic  methods  than  by  the 
beam- theory  and  strip-theory  methods  commonly  enqployed  In  the  past. 
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This  report  illustrates  the  development  and  application  of  a  method 
of  flutter  analysis  which  takes  into  account  three-dimensional  structural 
and  aerodynamic  behavior.  In  treating  the  structural  problem,  the  flutter 
mode  may  be  approximated  by  a  series  of  either  natural  or  assumed  vibra¬ 
tion  modes  which  could  have  platelike  distortions  or  shapes.  In  treating 
the  aerodynamic  problem,  lifting-surface  theory  is  used  to  obtain  aero¬ 
dynamic  forces  which  take  into  account  finite  span  and  cooqpressibility 
as  well  as  the  modes  of  vibration  of  the  structure. 

In  order  to  illustrate  the  application  of  the  method,  a  number  of 
flutter  calculations  sure  performed.  Primary  attention  is  directed 
toward  correlation  of  the  calculated  result  with  an  experimental  flutter 
result  for  a  delta  semi  span  wing  with  a  leading-edge  sweep  angle  of  ^5° 
at  a  Mach  number  of  0.8^.  The  calculations  are  bsised  on  carefully  meas¬ 
ured  natural -vibration  modes  obtained  by  means  of  an  optical  method.  A 
detailed  description  of  the  determination  of  these  modes  and  of  the  mass 
distribution  is  given  in  appendix  A. 

Other  ccLlculated  results  are  presented  for  the  same  wing  plan  form 
in  order  to  show  some  effects  of  variations  in  air  density  and  in  Mach 
number  (for  a  Mach  nvunber  range  from  0  to  0.95) •  Related  questions  con¬ 
cerning  the  number  of  modes  required  for  convergence  and  certain  effects 
of  nonorthogonality  of  the  measured  modes  are  considered.  Numerical 
evaluation  of  the  elements  of  the  flutter  determinant  is  discussed  in 
appendix  B. 


SIMBOIS 


generalized  aerodynamic  force  (see  eq.  (8)) 


arbitrary  constant  in  series  form  of  pressure  dis¬ 
tribution  (see  eq.  (l?)) 


b  local  wing  semichord,  ft 

bQ  wing  root  semichord,  ft 


d  distance  between  mirrors  on  wing  siirface  and  screen 

for  optical  method  of  measuring  modes,  in. 

Fn(6)  function  based  on  chordwise  pressure  term  (see 

eq.  (B2)) 

f  circular  frequency,  cps 
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^nm 

g 

h(x,y,t) 


constants  used  In  numerical  Integration  (see 
table  VI) 

dlnensloaless  pressure  mode  (see  eqs.  (17)  and  (I8)) 

structural  dancing  coefficient 

Instantaneous  deflection  of  point  on  surface 
in  the  flutter  node,  ft 


k 


normalized  displacement  at  point  x,y  In  ith  mode  of 
vibration,  hi(x,y) 


Ic 


chordwlse  Integrating  factor  appropriate  to  station 
(see  eq.  (B5)) 


Is 


spanvlse  Integrating  factor  appropriate  to  station  s 
(rice  eq.  (b6)) 


I 


(i) 

nzn 


K(M,ko,x-^,y-q) 


surface  integral  in  generalized  aerodynamic  force 
(see  eq.  (B2)) 

kernel  function  of  three-dimensional  Integral  equa¬ 
tion  (see  eq.  (12)),  l/sq  ft 

dimensionless  form  of  three-dimensional  kernel  func¬ 
tion  (see  eq.  (ik)) 

reduced-frequency  parameter,  bom/V 


dimensionless  pressure  function  (see  eq.  (5)) 


I 

M 

Ml 


Mlj 


vlng  semispan,  ft 
Mach  number 

generalized  mabs  associated  vlth  1th  mode  of  vibra¬ 
tion  (see  eq.  (3))»  slugs 

generalized-mass  coupling  term  (see  eq.  (21)),  slugs 


i*(3c,y) 


loceQ.  mass  per  unit  area  at  point  x,y,  slugs/sq  ft 

nuBiber  of  sheets  of  aluminum  foil  on  top  or  bottom 
surface  of  vlng 


N 


h 
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^(x,y,t) 


q 


i 


s 

■t 

V 

Vv 

x^y,5,n 

X,? 


y.n 


z 

a. 

ex' 


ex, a' 


local  oscillating  pressure  difference  between  top 
and  bottom  svirfaces  of  wing  in  the  flutter 
mode  (see  eq.  (4)),  Ib/sq  ft 

pressure  difference  at  point  x,y  in  Jth  mode  of 
oscillation,  ^j(x,y) 

generalized  coordinate  in  1th  mode  of  oscillation, 


conqplex  anqplltude  of  generalized  coordinate  in 
1th  mode,  ft 

area  of  wing  surface,  sq  ft 
time,  sec 

thicknesses  of  aluminum  Insert,  :4lumlnum  foil,  and 
balsa,  respectively  (see  eq.  (Al)) 

velocity  of  alrstream,  fps 

weight  of  wing  per  unit  area,  Ib/sq  in. 

Cartesian  coordinates  (see  sketch  following  eq.  (l2)) 

dimensionless  chordwise  variables  referred  to  bQ 
(see  eq.  (l4)) 

dimensionless  spanwlse  variables  referred  to  I 
(see  eq.  (l4)) 

airfoil  ordinate  (see  eqs.  (A2)  and  (A5)) 

local  slope  of  wing  in  pitch  during  vibration,  dh/dx 

local  slope  of  wing  in  a  direction  noxmal  to  line 
(or  ray)  of  constant  percent  chord 

vectors  representing  angular  displacements  a  and 
a* ,  respectively,  according  to  the  rlg^t-hand 
vector  rule  (see  eqs.  (A5)  and  (a6)  and  fig.  7) 

unit  weights  of  aluminum  alloy  and  laminated  balsa, 
respectively  (see  eq.  (Al)),  Ib/cu  in. 

reflected  displacement  on  screen  in  optical  method 
for  measxiring  mode  shapes,  in. 
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streanwise  and  epanvise  cooponents,  respectively, 
of  6 


e 


P 

0 


0) 

Indices: 

C,6 

f 


n,m 


r 

X 


y 


angtilar  chordwise  variable  (see  eq.  (l9))>  deg 
angle  of  sveepback  for  rth  ray,  deg 

generalized  mass-density  ratio  in  ith  mode  of  vibra¬ 
tion  (see  eq.  (lO)) 

air  density,  slugs/cu  ft 

loced  slope  of  wing  in  roll  during  vibration,  bh/^ 

local  slope  of  wing  in  a  direction  parallel  to  line 
(or  ray)  of  constant  percent  chord 

vectors  representing  angiilar  displacements  and  0', 
respectively,  according  to  the  right-hand  vector  rule 

angle  of  Image,  referred  to  the  horizontal,  on  screen 
In  optical  method  for  measuring  mode  shapes,  deg 

conplex  eigenvalue  of  determinantal  flutter  equation, 

{#'  •  W 

angular  frequency,  2jif,  radians /sec 


chordwise  and  spanwise  stations,  respectively 

flutter 

modes  of  vibration  under  consideration 

chordwise  and  spanwise  pressiu’e  modes,  respectively, 
in  aerodynamic  quantities  (m  denotes  power  to 
which  spanwise  variable  is  raised;  see  eq.  (18)) 

line  (or  ray)  of  constant  percent  chord  on  wing 

streamwise  direction  (related  to  local  pitching 
direction) 

spanwise  direction  (related  to  local  rolling 
direction) 


6 


NACA  TO  U595 


Matrix  notations: 


rectangular  matrix 
row  matrix 

column  matrix 
diagonal  nxtrlx 

Dots  over  symbols  denote  differentiation  with  respect  to  time. 


METHOD  OF  FUfTHR  ..tlALJfSIS 


The  present  section  Is  concerned  with  the  development  of  a  working 
form  of  the  flutter  equations  of  a  Baylelgh-Rltz,  or  modal,  type  of 
analysis  based  on  aerodynamic  forces  obtained  by  subsonic  lifting-surface 
theory.  A  brief  discussion  of  the  method  of  obtaining  the  aerodynamic 
forces  is  Included. 


Development  of  Flutter  Equations 

In  the  development  of  the  flutter  equations,  a  basic  assusptlon  of 
the  Rayleigh-Ritz  approach  is  that  the  displacement  h(x,y,t)  corre¬ 
sponding  to  the  flutter  mode  may  be  represented  by  a  superposition  of 
either  natural  or  assumed  modal  functions  In  the  fom 


»  q^hi  +  q2h2  +  .  .  .  (l) 

where  qj^  *  qj[e^  Is  the  generalized  coordinate  of  the  ith  chosen  mode 
and  hi  »  hi(x,y)  Is  the  corresponding  normalized  mode  shape. 

In  the  present  Investigation  the  nonuzlized  modsO.  functions  hi 
SLre  chosen  as  the  natural  (orthogonal)  modes  of  vibration  of  the  struc- 
tiore  under  consideration.  With  this  choice  of  mode  shapes,  the  general 
equation  of  motion  (obtained,  for  example,  as  in  ref.  1)  In  the  1th 
degree  of  freedom  may  be  written  In  the  form 


Miqi  +  a)i%Iiqi  «  JJ hi  Z^(x,y,t)  dS 


S 


(2) 
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vhere  Mi  represents  the  generalized  mass  in  the  1th  mode  and  is  defined 
as 


Ml  » 


(3) 


The  Integral  on  the  right  side  of  equation  (2)  represents  the  generalized 
aerodynamic  force  and  contains  the  aerodynamic  loading  iJ^(x,y,t)  which, 
consistent  with  equation  (l),  is  regarded  as  a  superposition  of 
aerodynamic-loading  inodes 


^(x,y,t)  =  i!^j^(x,y)  +  qg  +  .  .  .  (4) 


where  ^j(x,y)  denotes  the  aerodynamic  loading  associated  with  the 

—  liut 

mode  shape  hj.  In  application,  with  q^  =  q^e  ,  ^(x,y,t)  Is 
expressed  in  terms  of  dimensionless  functions  Li  through  the  relation 


Mx,y,t) 


(5) 


The  form  of  the  conplex  functions  Lj  is  dealt  with  at  a  later  stage. 

At  this  stage,  use  may  be  made  of  equations  (l)  and  (5)  to  obtain  equa¬ 
tion  (2)  in  the  form 


Miqi 


4xpV^Z 

bo^ 


JJ  hiC^Li 


+  q2^  •  •  •  )dS  =  0  (6) 


S 


or,  alternatively. 


i-P 

\  CO/ 


^2^2 


.) 


Miko 


(7) 
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where 


and 


(8) 


If  equation  (7)  Is  written  for  each  of  the  degrees  of  freedom  imder 
consideration,  a  set  of  sisultaneous  honogeneovis  equations  results.  The 
flutter  condition  is  then  given  by  the  vanishing  of  the  determinant  of 
the  coefficients  of  ^  in  these  equations;  thus, 


‘11 


1-i— « 

A21 


‘12 


1  - 


—  n  + 


^22 


(9) 


where 


1 


(10) 


and  where  A  is  defined  by 


n 


(1  +  ig) 


(11) 


In  this  form  (1  is  a  c(Mq>lex  eigenvalue  and  contains  the  unknown  fre¬ 
quency  (u  and  a  dB]iq>lng  coefficient  g  that  becones  tero  at  the  border¬ 
line  flutter  condition  between  daaqped  and  undaiiq>ed  motion. 
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Equation  (9),  together  with  the  definitions  given  in  equations  (3)^ 
(8),  (10),  and  ^l),  constitutes  the  working  form  of  the  flutter  equa¬ 
tion.  In  other  sections  of  the  report  the  determination  of  the  Ingre¬ 
dients  of  eq\iation  (9)  is  discussed;  appendix  A  Is  concerned  with  the 
determination  of  mode  shapes  and  mass  distribution  for  a  particular  con¬ 
figuration,  and  e^pendlx  B  deals  with  numerical  techrdques  used  to 
evaluate  generalized  masses  and  generalized  aerodynamic  forces. 


Determination  of  Aerodynamic  Loading 

In  the  determination  of  the  aerodynamic  loading,  the  functions  Lj 

required  in  the  generalized  aerodynamic  forces  are  obtained  from  the 
integral  equation  which  relates  lift  and  downwash  distributions  in  sub¬ 
sonic  lifting-stirface  theory.  (See,  for  eican5)le,  ref.  2.)  A  systematic 
numerical  solution  is  employed  herein  which  is  similar  to  that  of  refer¬ 
ence  3  but  which  makes  use  of  a  more  exact  form  of  the  kernel  function 
together  with  more  refined  miner Ical-integratlon  procedures.  Qhe  method 
enqployed  has  been  programed  for  the  IBM  type  704  electronic  data  pro¬ 
cessing  machine. 

For  the  purpose  of  the  present  investigation,  the  Integral  equa¬ 
tion  may  be  written  as 


*(l^  *  *  5;^  «  dl)  (12) 


s 


where  the  coordinate -axis  system  is  shown  in  the  followlzig  sketch: 
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where  and  denote,  respectively,  the  coordinates  of  the 

leading  and  trailing  edges.  The  term  on  the  left  side  of  equation  (12) 
denotes  the  downwash  angle  at  point  x,y  associated  with  a  dlsplsicenent 
h(x,y,t).  On  the  right  side  of  eqxiatlon  (12),  the  kernel  function 

K^M,~,x-|,y>qj  represents  the  downwash  produced  at  point  x,y  by  a  unit 

pressure  load  at  i>olnt  |,t);  the  function  ^(|,q,t)  represents  the  load 
required  to  satisfy  the  specified  downwash  condition.  (!nie  minus  sign 
on  the  left  side  of  eq.  (12)  Is  associated  with  the  fom  of  the  kernel 
function  and  4rlses  from  the  use  herein  of  the  sign  convention  for  down- 
wash,  displacement,  and  lift  which  are  positive  downward.) 

Substituting  equations  (l)  and  (3)  Into  equation  (l2)  gives 


♦ 


Introducing  the  following  dimensionless  quantities 


7  =  ^ 

,.a 


K(M,ko,x-r,y-n) 


l2  K| 


^M,tt,x-|,y-T\^ 


(14) 
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Into  equation  (ij)  gives 


From  equation  (15)  the  form  of  the  integral  equation  for  the  Jth 
xnode  of  osclUatica  may  be  seen  to  be 


x-Ly-n) 


Lj  d|  dn 


(16) 


In  the  solution  of  equation  (16),  it  Is  assumed  that  a  function  Lj  may 
he  represented  by  a  series  of  pressure  modes  of  the  form 

Lj  «  ^  ^  ^nm 
n  m 

vhere  the  quantities  are  aibitrary  constants  to  be  determined. 

The  indices  n  and  m  are  associated  with  chordvlse  and  spanwise  pres¬ 
sure  modes,  respectively.  The  form  of  the  functions  fmn  Is  dictated 

by  known  leading-  and  tr ailing-edge  conditions.  In  subsonic  flow  they 
sho\ild  satisfy  the  edge  conditions  pertinent  to  the  Kutta  condition.  In 
the  existing  procedure  for  solving  the  integral  equation  they  have  been 
expressed  (in  terms  of  an  angular  chordwlse  variable  0)  as 
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vhere,  for  exai^le,  for  synnetrlc  motion  the  Index  m  Is  equal  to 
0,  2|  4,  .  .  .  and  where  variable  6  (see  sketch  vhlch  follows 
eq.  (12))  Is  related  to  i  by  the  expression 

I  .  .11.^03  8  (19) 
2  bo 


in  which 


^te  ^le 
2 


denotes  the  equation  of  the  wing  midchord. 


If  the  series  form  of  Lj  given  by  equation  (17)  is  substituted 

Into  equation  (l6),  a  working  form  of  the  Integral  equation  Is  obtained 
as 


JJ'fnm  K(M,ko,x-r,7-n)  d?  dn 


8 


(20) 


and  may  be  seen  to  consist  of  a  sxuanatlon  of  definite  surface  integrals, 
each  weighted  by  an  unknown  constant  a^^ . 

In  order  to  detexnlne  the  constants  a^,  a  collocation  procedure 

Is  used.  The  right  side  of  equation  (20)  Is  evalviated  for  as  many 
points  x,y,  designated  as  control  points,  as  there  are  unknown 
constants  Onm  1*  equated  to  the  known  downwash  angle  at  each  con¬ 
trol  point.  A  set  of  simultaneous  equations  Is  thus  obtained  which  may 
be  solved  for  the  values  of  aQ^.  Once  determined,  the  constants  ag^ 
are  used  with  equations  (l7)  and  (l8)  to  define  the  pressure  functlcm  Lj; 
and  Lj  may  then  be  en^loyed  In  equation  (8)  to  obtadn  the  generalized 
aerodynamic  forces  Ajj.  In  the  present  a^llcatlon  nine  terms  of  the 

series  In  equation  (17)  were  used  (so  that  n  «  0,  1,  2  and  m  >  0,  2,  4), 
and  the  downwash  was  satisfied  at  the  nine  control  points  shown  In  fig¬ 
ure  l(a).  As  previously  noted,  a  solution  of  the  Integral  equation  as 
represented  by  equation  (20)  has  recently  been  programed  for  the  OH  type 
704  electronic  data  processing  machine. 


APPUCATION  OF  METHOD  TO  A  DELTA  WSfO 


In  order  to  show  the  application  of  the  analytical  treatment  described 
in  the  preceding  section,  a  number  of  calculations  have  been  made  for  a 
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delta  semlspan  vlng  %rlth  a  leading-edge  sweep  angle  of  Prljaary 

a'ttentlon  Is  directed  toward  a  correlation  of  calculated  results  with 
an  experlnental-flutter  result.  Other  calculations  are  perfoxned  for 
the  configuration  of  the  experiment  to  explore  some  effects  of  varying 
Mach  number  and  air  density.  In  addition^  ^one  effects  of  vazying  maas 
eund  natural  frequencies  and  certain  effects  of  nonorthogonality  of  the 
neasiired  nodes  are  examined. 


Correlation  of  Calculations  With  Experiment 

The  experimental  result  used  as  a  basis  for  the  calculations  has 
been  obtained  by  William  T.  Lauten,  Jr.^  and  Marvin  F.  Burgess  at  the 
Langley  Aeronautical  Laboratory.  Construction  details  of  the  delta 
semlspan  wing  under  consideration  are  shown  in  figure  2.  The  mauts  prop¬ 
erties  of  the  model  are  shown  as  chordwlse  and  spanwlse  distributions 
of  weight  per  vmlt  area  In  figure  3*  The  mode  shapes  In  the  first  foxor 
natural  modes  of  vibration  were  obtained  by  means  of  an  optlcsil  method 
and  sure  shown  In  figure  4.  Detedls  of  the  methods  used  In  obtaining 
mass  and  node  shapes  are  described  In  appendix  A. 

Tvo  models  were  Involved  In  the  experimental  program.  One,  des¬ 
ignated  as  model  A,  fluttered  to  destruction  before  Its  mode  shapes  and 
mass  distribution  had  been  determined.  These  structural  properties 
were  obtained  by  use  of  a  second  model,  model  B,  which  was  built  to  the 
sane  specifications  and  had  very  similar  nodal  patterns.  As  can  be  seen 
In  table  I,  model  B  was  lighter  than  model  A  and  had  different  natural 
frequencies.  In  all  flutter  calculations  the  node  shapes  for  model  B 
were  used,  and  In  the  calculations  for  model  A  the  generalized  masses 
were  adjusted  by  the  ratio  of  the  total  masses  of  the  models. 

The  flutter  calculation  of  primary  Interest  has  been  made  for 
model  A,  vhlch  fluttered  as  a  cantilever  at  a  Mach  nuniber  of  O.63  with 
eui  air  density  of  O.OOO787  slugs/cu  ft.  Results  of  this  flutter  cal¬ 
culation  are  shown  in  table  II(a}.  Converged  flutter-speed  solutions 
sure  obtained  (as  Indicated  by  the  agreement  of  the  first  three-  or  four- 
mode  calculations)  which  are  about  3  percent  less  than  the  experimental 
walue.  Calculated  results  based  on  two  twonnode  subcases  are  Included 
In  table  Il(a).  The  flutter  speed  calculated  with  modes  1  and  3  Is 
about  2  percent  above  the  experimental  flutter  speed.  It  Is  noted  that 
modes  1  and  3  resemble,  respectively,  flrst-bendlng  and  first-tors  Ion 
modes  of  a  beam. 


Variations  In  Mach  Number  and  Air  Density 

Results  of  converged  flutter  calculations  at  other  values  of  Mach 
nunber  and  air  density  for  the  configuration  of  the  experimental  flutter 
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condition  are  given  In  table  II (b)  and  are  presented,  In  different  forms. 

In  figures  3  and  6.  Figure  3  Illustrates  the  effect  on  Vf  of  varying 

density,  with  Mach  number  held  constant  at  0.83>  Results  are  shown  In 
terms  of  an  alr-denslty  peurameter  l/^p.  The  effect  of  Mach  number  on 
calculated  f.luttf>r  speed  Is  shown  In  figure  6  for  each  of  three  values 
of  air  densl  Results  are  presented  In  terms  of  a  flutter-speed 

coefficient  \f/b(ya^.  As  can  be  seen,  the  flutter  speeds  calculated  for 

the  present  configuration  appear  to  remain  nearly  constant  for  all  Mach 
numbers  up  to  M  *  0.93*  The  corresponding  calculated  flutter  frequencies, 
on  the  other  hand,  tend  to  decreaise  with  Increasing  Mach  number,  as  may 
be  seen  In  table  Il(b). 


Effects  of  Certain  Structural  Modifications 

The  effects  of  the  differences  In  total  mass  euid  natured.  frequencies 
between  models  A  and  B,  noted  In  teible  I,  were  considered  worthy  of  fur¬ 
ther  study.  First,  flutter  cflilculatlons  were  performed  for  model  B,  and 
the  converged  flutter  speed  was  10  percent  less  than  that  for  model  A. 

(Ccmpare  case  1  of  table  Il(c)  with  the  four^node  result  In  table  11(a).) 

In  an  effort  to  separate  the  effect  of  the  difference  In  total  mass  from 
the  effect  of  the  difference  In  natural  frequencies,  additional  flutter  ,  ^ 

calculations  based  on  the  generalized  masses  of  model  B  and  the  natural 
frequencies  of  model  A  were  performed.  The  resxilts,  listed  In  table  Il(c) 

(case  2),  show  that  about  three-fourths  of  this  10-percent  (Ufference  In  .  < 

flutter  speed  was  due  to  the  differences  In  the  natural  frequencies 
between  the  two  models.  The  remaining  difference  In  flutter  speed  Is 
then  attributed  to  the  difference  in  total  mass. 

As  can  be  seen  In  equation  (9)/  the  off diagonal  elements  of  the 
flutter  determinant  contain  no  inertial-coupling  terms.  This  is  a  con¬ 
sequence  of  the  orthogonality  condition  for  natural  coupled  modes,  which 
may  be  expressed  as 


MlJ  » 


(i  /  J)  (21) 


It  is  recalled  that,  in  the  present  analysis,  measured  modes  were 
employed.  In  order  to  Investigate  their  orthogonality,  values  of  M^j 
were  conqputed  for  model  A  by  use  of  eqviation  (21)  and  are  listed  in 
table  III.  In  order  to  determine  the  effects  of  the  nonzero  values 
of  Mij,  they  were  added  to  the  appropriate  elements  of  equation  (9) 
and  two  cal cvilat ions  were  made  for  model  A.  First,  the  aerodynamic 
terms  were  eliminated  and  the  natiiral  frequencies  were  cong>utedj  results 
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are  shown  In  table  III  and,  with  the  exception  of  are  seen  to  differ 

by  a  negligible  amount  frcxn  those  in  table  I.  In  a  second  calculation  the 
change  in  flutter  speed  due  to  including  the  terms  was  determined; 

the  other  parameters  used  were  those  of  the  four-mode  solution  of 
table  11(a).  Results  are  listed  as  case  3  of  table  Il(c)  and  show  that 
the  converged  flutter  speed  has  been  decreased  by  less  than  1  percent. 


CONCLUDING  REllARKS 


The  development  and  the  numerical  application  have  been  presented 
of  a  method  of  flutter  analysis  which  takes  into  account  three-dimensional 
structural  and  aerodynamic  behavior.  The  flutter  mode  was  approximated 
by  a  series  of  natxiral -vibration  modes,  and  the  aerodynamic  forces  cor¬ 
responding  to  these  modes  were  derived  from  subsonic  lifting-surface 
theory,  according  to  the  kernel-function  approach,  for  a  finite  wing 
oscillating  in  con5)ressible  flow. 

The  application  was  made  to  a  delta  semlspeui  wlxig  with  a  leading- 
edge  sweep  angle  of  45°  which  fluttered  at  a  Mach  number  of  O.85. 

Results  of  flutter  calculations  show  that  \rtien  the  first  three  or  four 
natural- vibration  modes  were  Included,  conver^d  solutions  for  the  flutter 
speed  were  obtained  which  were  about  5  percent  less  than  the  experimental 
value.  Further  application  of  the  analysis  was  made  to  orient  the  experi¬ 
mental  flutter  condition  with  theoretical  flutter-speed  boundaries  for  a 
range  of  densities  suid  Mach  numbers  including  those  of  the  e;q>eriment  and 
to  study  the  effects  of  certain  variations  in  natural  frequencies  and 
total  mass. 


Langley  Aeronautical  Laboratory, 

National  Advisory  Committee  for  Aeronautics, 
Langley  Field,  Va.,  July  I7,  1958 • 
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APPENDIX  A 

STRUCIURAL  PROPERTIES  OF  DELTA-WING  MODEL 


The  structural  properties  of  the  delta-ving  model  needed  for  the 
flutter  analysis  that  is  presented  and  applied  in  this  report  consist 
of  the  mass  distribution  and  the  natural  modes  of  vibration,  vdiich  are 
represented  in  equations  (3)  and  (8)  as  m(x,y)  and  respectively. 

This  appendix  provides  a  description  of  the  methods  used  to  determine 
these  properties  for  the  delta-ving  configuration  considered  in  this 
analysis. 


Description  of  Model 

Figure  2  shows  pertinent  details  of  the  construction  of  the  model. 
This  construction  consisted  of  a  0.02-inch- thick  sheet  of  aluminum  alloy 
located  at  the  plane  of  symmetry  vith  vertically  laminated  balsa  glued 
to  each  surface  of  the  sheet,  nie  balsa  laminations  were  shaped  to  the 
ordinates  of  an  NACA  65A004  airfoil  section,  and  the  outside  surface  was 
wrapped  vith  layers  of  0.001-lnch-thlck  aluminum  foil.  The  wing  tip  was 
rounded  so  that  the  semi  span  was  32.63  Inches  at  a  point  1.3  inches  for¬ 
ward  of  the  trailing  edge.  The  aerodynamic  aspect  ratio,  with  the  assunq)- 
tion  of  a  straight  or  squared-off  tip  based  on  extensions  of  leading  and 
trailing  edges  to  a  common  point,  was  3*34. 

Two  models  were  Involved,  both  built  to  the  same  specifications. 

One  model,  designated  as  model  A,  fluttered  to  destruction  before  all  of 
its  properties  were  determined;  the  other  model,  model  B,  was  used  to 
obtain  the  mass  distribution  and  experimental  natural^node  shapes  for 
these  elements.  Certain  mechanical  properties  of  both  models  are  com¬ 
pared  in  table  I. 


Mass  Properties  of  Model 

The  mass  properties  of  both  models  studied  in  this  investigation 
are  based  on  the  chordvlse  and  spanvlse  distributions  of  weight  per  unit 
area  shown  in  figure  3*  These  distributions  were  calculated  from  the 
equation 

''v' ’’bS 

where  v^  is  the  weie^t  of  the  wing  per  unit  area,  7^  is  the  unit 
weight  of  aluminum  alloy,  tj  is  the  thickness  of  the  aluminum  Insert, 
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tp  is  the  thickness  of  the  alumlnun  fell,  tfi  Is  the  thickness  of  the 
laainated  balsa,  emd  is  the  unit  vei^t  of  laminated  balsa.  Ihe 
value  of  7q,  \diich  was  determined  from  the  measured  total  veiedit  of  the 

ving  and  the  volvune  occupied  by  the  balsa,  was  foimd  to  be  0.01CX)7  Ib/cu  in. 
By  using  this  value  together  with  y^  *  0.100  Ib/cu  in.  and  the  data  given 

in  figure  2,  equation  (Al)  can  be  written  in  the  form 

«  0.100(0.020  +  0.002N)  +  0.01007[48zb  -  (0.02  +  0.002N)j  (A2) 

everywhere  on  the  wing  except  in  the  regions  of  the  leading  and  trailing 
edges  which  were  covered  by  an  additional  layer  of  aluminum  foil  which 

was  inches  wide  on  each  surface.  Hie  consequent  increase  in  weight 
in  these  edge  regions  was  accounted  for  in  the  equation 

Wy  »  0.100(0.022  +  0.002N)  +  0.01007[48£b  -  (0.022  +  0.002N)]  (A3) 

In  equations  (A2)  and  (A3)  N  is  the  number  of  0.001-lnch-thick 
sheets  of  aluminum  foil  on  each  (top  or  bottom)  surface  of  the  wing,  and 
z  is  the  ratio  of  one-half  the  local  thickness  to  the  local  chord  for 

an  NACA  63A004  airfoil. 


Optical  Method  of  Measuring  Natursd-Vibration  Modes 

The  natured-vlbration  modes  were  obtained  from  the  results  of  shake 
tests  conducted  at  zero  airspeed  and  involving  the  use  of  an  electro- 
dynamic  shaker.  The  shaker  was  located  near  the  wing  root  in  order  to 
reduce  as  much  as  possible  the  effect  of  the  shaker  mass  on  the  natural 
modes.  For  each  natural  frequency  found,  the  corresponding  mode  shape 
up  throuf^  the  fourth  mode  was  determined  by  application  of  an  optical 
method  described  in  this  section.  This  method  Involved  the  measurements 
and  numerical  integration  of  local  slopes  in  pitch  and  roll  directions 
at  24  stations  distributed  over  the  wing  surface. 

A  schematic  'diagram  of  the  apparatus  used  in  this  optical  method  is 
shown  in  figure  J.  Small  mirrors  that  were  l/l6  inch  square  were  glued 
to  the  wing  at  six  stations  along  rays  at  23,  50,  and  75  percent  chord 
and  along  the  trailing  edge  as  shown  in  the  mirror  lattice  in  the  upper 
left  corner  of  the  figure.  Light  aimed  at  the  wing  was  reflected  by  the 
mirrors  onto  a  screen  located  a  considerable  distance  from  the  model. 
During  vibration  in  a  natural  mode.  Images  reflected  from  these  mirrors 
appeared  as  straight  lines  whose  lengths  euid  directions  were  marked  on 
the  screen.  “Riese  measured  lengths  and  directions  were  resolved  into 
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local  pitching  and  rcUlng  slopes  for  small  values  of  a  and  ^  by 
the  relations 


(Ak) 


It  should  be  noted  that,  although  this  relation  is  shown  for  an  instan¬ 
taneous  position  of  the  model  from  equilibrium,  it  also  applies  when  6 
corresponds  to  a  double  amplitude,  as  is  the  case  during  vibration.  The 
distance  d  in  this  application  was  155  inches. 


Ihe  local  slopes  in  pitch  and  roll  were  obtained  at  six  stations 
along  each  of  four  rays  at  25,  50,  75,  and  100  percent  chord.  Deflec¬ 
tions  along  each  ray  were  determined  by  numerically  Integrating  spanwise 
curves  faired  throu^  the  six  values  of  along  the  ray.  The  vec¬ 

tors  a,  a',  and  indicated  in  the  lower  left  comer  of  fig¬ 
ure  7  represent  angular  displacmnents  a,  a',  0,  and  0',  respectively, 

according  to  the  right-hand  vector  rule,  and  from  this  vector  relationship 


a  =  a*  cos  Aj.  +  0*  sin  Aj. 
0  =  -a'  sin  Aj.  +  0'  cos  Aj. 

from  which 


(A5) 


0'  =0  cos  Aj,  +  a  sin  Aj. 


(A6) 


Thus,  for  the  ith  natural  mode,  after  substitution  of  equation  (A4)  into 
equation  (a6). 


^  *  I'S.S,'*)  (A7) 

where  the  subscripts  s  and  r  identify  a  point  on  the  mirror  lattice 
as  indicated  in  the  upper  left  corner  of  figure  7.  Values  of  0'^^^ 

8Z* 

determined  from  this  equation  are  given  in  table  IV.  For  modes  higher 
than  the  first  mode,  the  sign  of  the  angle  ^4-  a^,  was  governed  by 

the  observed  node  locations  on  the  rays  rather  than  by  the  actiial  numer¬ 
ical  value  of  this  angle. 
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iniese  slopes  were  in  t\irn  used  to  determine  the  deflections  by 
numerically  integrating  spanwise  curves  faired  through  the  six  values 
of  along  each  ray.  A  check  of  the  data  disclosed  that,  for  any 
point  on  a  particular  ray,  the  deflection  obtained  In  this  manner  was 
in  very  close  agreement  with  the  deflection  obtained  by  integrating  0' 
along  the  trailing  edge  to  the  spanwise  position  of  the  point  and  then 

by  integrating  the  curve  faired  through  four  chordwise  values  of 

forward  of  the  trailing  edge  to  the  i>oint  itself.  On  the  basis  of  this 
check,  the  normalized  deflection  curves  shown  in  figure  4  for  the  first 
four  modes  have  been  adjusted  to  be  coopatlble  with  the  measured  values 

of  given  in  table  V. 

The  integration  of  along  a  ray  was  performed  by  means  of  the 
direct- sumnatlon  method  Illustrated  in  reference  4  with  the  use  of  10 
equally  spaced  stations  for  the  first  mode,  15  to  17  stations  for  the 
second  mode,  11  to  l4  stations  for  the  third  mode,  and  l8  and  19  sta¬ 
tions  for  the  fourth  mode. 

•Hie  nodes  shown  in  figures  4(b),  4(c),  and  4(d)  trace  out  very 
nearly  the  same  paths  observed  on  model  A,  particularly  in  the  outboard 
region  of  the  wing  plan  form.  Because  this  similarity  is  a  good,  though 
not  cookie te,  indication  of  mode-shape  similarity  between  the  two  models, 
the  mode  shapes  used  to  calculate  the  generalized  masses  of  model  A  were 
assumed  to  be  the  same  as  those  measured  for  model  B. 
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APPENDIX  B 

EVALUATION  OF  EliMENTS  OF  DETENHIIUNTAL  FLUTTER  EQUATION 
AND  APPLICATION  TO  A  SPECIFIC  CONFIGURATION 


Evaluation  of  Generalized  Aerodynamic  Force 

In  order  to  evaluate  the  generalized  aerodynamic  force  Aj^j 

(eq.  (8)),  use  Is  made  of  the  dimensionless  variables  introduced  by  equa¬ 
tion  (l4)  and  the  pressure  function  Lj  given  by  equations  (17)  to  (I9) 

to  write  equation  (8)  as 


*  folo  l(' 


^00  ^  ^02 


f 


sin  e (a^o  +  ?  a-j^  +  .  .  .  )  + 


12 


^  sin  26^1 


ago  +  y%2 


■) 

.  .  •  •  'Jsln  0  d0  dy 


(Bl) 


In  the  evaluation  of  this  equation,  the  following  definitions  are  useful: 


Fq(0)  =  cot  ^  sin  0=1+  cos  0 

F;,^(0}  =  sln^0 

F2(0)  =  ^  sin  20  sin  0 

\ 


(B2) 


Equation  (Bl)  may  be  written  In  terms  of  the  Integral  I^^  as  the 
summation 


IZ 

n  m 


(B3) 


HKCk  *ni  4395 


21 


or,  as  the  matrix  product 


•  •  -I 

ij  L  ™® J  I  J  L  J 


Numerical  integrating  techniques  can  be  applied  to  evaluate  the 
surface  Integrals  Ijy^  in  the  row  matrix  of  equation  (b4).  ^It  is 

noted  in  this  connection  that  the  use  of  the  angular  coordinate  0  in 
the  chordwise  integration  is  particularly  convenient  for  this  applica¬ 
tion.  For  exanqple,  when  represented  in  Cartesian  coordinates,  the  first 

term  in  the  integrand  of  is  infinite  at  the  wing  leading  edge: 

a 

when  expressed  in  polar  coordinates,  the  product  cot  ^  sin  0  is  a  non- 

singular  function.  )  A  number  of  spanwise  stations  may  be  chosen  at 
y  =  y^  with  appropriate  spanwise  integrating  factors  Ig.  Similarly, 
at  each  spanwise  station,  a  number  of  chordwise  stations  are  selected  at 
e  *  0g  with  appropriate  chordwise  integrating  factors  I^j  s  and  c 
thus ' identify  stations  on  an  integrating  lattice  such  as  that  shown  in 
figure  1(b) . 


A  matrix  product  of  chordwise  terms  may  be  formed  as 


^o(»2)  • 

• 


and  a  matrix  product  of  spanwise  terms  may  be  formed  as 
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The  normalized  displacements  at  the  integrating  stations  may  be 

arranged  in  a  rectangular  matrix  of  the  form 


^i(®l*^l)  ^i{®l"y2)  *  •  •  ^i(®l^ys) 

hi(02,yj 


(B7) 


||^i(®c^yi) 

Equations  (B5)  to  (B7)  nay  be  used  to  obtain  the  surface  integrals 

/  \ 

of  equation  (B2)  from  the  matrix  product 

m-* 

=  [«][»l]PJ  (B8) 


The  elements  of  the  rectangular  matrix  on  the  left  side  of  equation  (b8) 
may  then  be  rearranged  as  the  row  matrix  equation  (b4)  and  be 

used  with  the  values  of  the  constants  a^*^^  to  obtain  the  generalized 


nm 


aerodynamic  force  Aj^j.  Equation  (b4)  may,  of  course,  be  expanded  for 
any  number  of  modes  into 


'll 


'21 


'12 


'll 


'ij 


00 


00 


nm 


00 


a 


(1) 

nm 


a(^) 

*00 


(B9) 
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Evaluation  of  Generalized  Mass 

Ihe  generalized  mass  Is  expressed  as  a  surface  Integral  by 

equation  (3)  and  may  be  evaluated  by  use  of  the  same  matrices  of  Inte¬ 
grating  factors  and  displacements  employed  In  evaluating  the  generalized 
aerodynamic  forces  A^j.  For  this  purpose,  equation  (3)  may  be  written 

in  terms  of  the  variables  defined  by  equations  (l4)  and  (19)  as 


1 1  rtli 


“  ^0^  f  J  ^1^  ^  sin  0  de  dy  (BIO) 

or,  as  a  sumnatlon  which  Involves  the  Integrating  factors  Ig  and  1^, 
used  with  Aj^j,  In  the  form 

Ml  “  V  X  Z  "(9c.ys)»i^(6e.yB)«ln  9^ 


s  c 

or  Ir  'Si.'ms  of  a  matrix  product  as 

c 


— —  ^  ^  f  ' 

s  i 


(B12) 


The  colximn  matrix 


M- 


which  pertains  to  a  particular  spanwlse  sta- 


bion,  may  be  obtained  by  sqxiaring  the  elements  in  the  corresponding 
colxnnn  In  the  matrix  J^hjJ  (eq.  (BT))*  ®ie  Inertial-coupling  term 

(eq.  (21))  is  found  by  replacing  the  matrix 


in  equa¬ 


tion  (B12)  with  a  matrix  of  crossproducts  ihj)  1  J  ■ 


Application  of  Method  to  a  Specific  Configuration 

In  the  application  to  a  specific  configuration  of  the  integrating 
procedxires  discussed  in  the  previous  section,  either  equal  or  unequal 
intervals  in  a  given  variable  may  be  taken.  With  equal  Intervals, 
appropriate  integrating  factors  are  readily  available j  when  unequal 
Intervals  are  employed,  corresponding  integrating  factors  must  be  devel¬ 
oped.  Useful  procedures  relating  to  numerical  integration  are  presented 
in  reference  5* 
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In  the  present  application  to  the  flutter  of  a  delta  wing,  10  equal 
Intervals  in  the  chordwise  variable  0  and  8  equal  Intervals  in  the 
spanwise  variable  y  were  used.  The  integration  stations  are  shcvn  in 
figure  1(b) .  Chordwise  integrating  factors  I^,  employed  are  based  on  a 

rule  employing  overlapping  quintic  functions  derived  in  reference  4. 
Spanwise  integrating  factors  Ig  were  based  on  a  rule  which  employs 

overlapping  quartlc  functions.  Data  pertinent  to  the  evaluation  of  the 
factors  Ij,  and  Ig  are  given  in  table  VI. 


Normalized  mode  shapes  hj^  employed  in  obtaining  both  generalized 
force  Aj^j  and  generalized  mass  were  obtained  from  the  plots  of 

figure  4.  TOe  normalizing  station  for  each  mode  was  arbitrarily  located 
on  the  midchord  at  y  »  0.875*  Values  of  mass  per  unit  area 

for  use  in  determining 

shown  in  figure  3* 


M4 


were  obtained  frcxi  the  weight  distribution 

TOe  ccsnponents  ^4-  and  h^  of  the  downwash  angle 

employed  as  a  boundary  condition  in  equation  (20)  are  listed  in  table  VII 
for  the  control  points  shown  in  figure  1(a) . 
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TABI£  I 

CERTAIN  MECHANICAL  PROPERTIES  OF  DELTA-WING  MODELS 


Item 

Model  A 

Model  B 

Aspect  ratio  (squared-off  tip)  . 

3.54 

3.54 

Airfoil  section  . 

NACA  65A004 

NACA  65A004 

Ay,  deg; 

Leading  edge  . 

45 

45 

25  percent  chord  (r=l)  . 

36.83 

36.83 

50  percent  chord  (r=2)  . 

26.58 

26.58 

75  percent  chord  (r=3)  . 

14.03 

14.03 

Trailing  edge  (r=4)  . 

0 

0 

^0^  . 

1.458 

1.458 

Z ,  f  t . 

2.739 

2.739 

Total  mass  of  wing,  Ib-sec^/ft  . 

0.1789 

0.1644 

Natural  frequencies,  cps; 

First  mode  . 

21 

21 

Second  mode  . 

58 

52.5 

Third  mode  . 

81 

77 

Fourth  mode  . 

115' 

107 

KACA  TN  4395 
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IAB1£  11 


RESULTS  OF  FLUTTER  ANALTSIS  AND  CORREIATICNI  WITH  EXPBRlNEirr 

(a)  Correlation  of  theory  and  experlaent  for  aodel  A 
(M  «  0.85;  P  •  0.000787  slugs /cu  ft) 


Nodes  used 
in  analysis 

Vf> 

fps 

Cp8 

^,t 

Results  of  calculations  . 

1,  3 

1,  2,  3 

1,  2,  3, 

2,405 

943 

882 

876.5 

73.5 

41.5 

40.0 

39.8 

0.280 

.403 

.415 

.416 

924 

37.9 

0.375 

(b)  Results  of  converged  flutter  calculations  at  other  Mach  numbers  and  densities 

(Four  modes  used  in  analysis) 


M 

slugs/cu  ft 

H 

m 

B 

slugs/cu  ft 

H 

^f' 

cps 

B 

0.85 

0.000504 

.001267 

.0021 

.00326 

39.3 

40.2 

42.5 

43.9 

0.329 

.511 

.630 

.726 

0 

.4 

.7 

.95 

0.000787 

916 

930 

920 

926 

53.5 

52.1 

46.0 

36.0 

0.535 

.513 

.458 

.375 

0 

.k 
•  7 
.95 

0.000504 

1,129 

1,133 

1,135 

1,073 

53.0 

51.2 

45.5 

36.9 

0.430 

.414 

.367 

.315 

0 

.4 

.7 

.95 

0.001267 

735 

741 

752 

792 

54.0 

52.6 

46.8 

42.2 

0.673 

.650 

.570 

.488 

(c)  Effects  of  structural  modifications 
(M  «  0.85;  P  =  0.000787  slugs/cu  ft) 


Case 

Model 

**1J 

(1  ftj) 

Results  of  converged 
solutions 

m 

m 

1 

B|B 

0 

797 

38.2 

0.439 

2 

0 

857 

40.0 

.426 

3 

B 

Honsero 

872 

39.0 

.410 

*The  model  combined  the  generalised  masses  of  model  B  and 
the  natural  frequencies  of  model  A. 


28 


NACA  TN  ^395 


TABLE  III 

GENERALIZED  MASSES  AND  CORRESPONDING  CALCUUTED  NATURAL 
FREQUENCIES  COMPUTED  FROM  EQUATION  (21)  FOR  MODEL  A 


■ 

1 

\y  Slugs,  for 

values  of  J 

of  - 

ops 

1 

2 

3 

4 

B 

0.010376 

-0.0015128 

0.0054925 

-0.026760 

20.95 

B 

-.0015128 

.0159856 

-.0139947 

-.0645656 

57.2 

B 

.0034925 

-.0139947 

1.245043 

-.0471072 

81.6 

B 

-.026760 

-.0645656 

-.0471072 

2.897119 

125.7 

NACA  TN  4395 
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TABIil  IV 


MEASURED  SLOPES  OF  NATURAL-VIBRATION  MODE  SHAPES 

sr 

CALCUIATED  BY  MEANS  OF  EQUATION  (A?)^ 


B 

s 

Semi  span 
station, 
percent 

r  =  1 

(25  percent 
chord) 

r  =  2 

(50  percent 
chord) 

r  =  3 

(75  percent 
chord) 

r  =  4 
(trailing 
edge) 

■ 

1 

20 

0.280 

0.727 

1.096 

1.350 

■ 

2 

35 

1.048 

1.570 

1.870 

2.550 

II 

5 

50 

2.298 

2.964 

3.410 

4.080 

B 

4 

65 

4.001 

4.495 

4.927 

5-284 

B 

5 

80 

5.842 

6.373 

6.755 

6.794 

■ 

6 

90 

7.707 

7.937 

7.916 

7.621 

1 

20 

-0.658 

-0.736 

-1.243 

-1.873 

2 

35 

-1.150 

-1.208 

-1.526 

-1.777 

5 

59 

-.993 

-.605 

-.175 

.481 

c. 

4 

65 

.992 

1.888 

2.868 

4.012 

5 

80 

6.029 

7.245 

8.683 

10.046 

6 

90 

11.938 

12.558 

12.978 

12.915 

1 

20 

0.891 

0.659 

-0.968 

-5.355 

2 

35 

1.037 

.065 

-1.262 

-3.117 

3 

50 

.783 

-.205 

-1.055 

-1.207 

5 

4 

65 

-.1195 

-.259 

-.116 

1.217 

5 

80 

-.622 

.125 

1.539 

3.636 

6 

90 

-.4505 

1.012 

2.854 

4.987 

■ 

H 

20 

0.436 

0.455 

0.667 

1.558 

■ 

35 

.303 

.239 

.243 

-.356 

B 

50 

-.591 

-.987 

-1.492 

-5.022 

B 

65 

-1.830 

-1.720 

-1.877 

-2.200 

B 

El 

80 

.373 

1.486 

5.017 

4.788 

■ 

90 

7.963 

9.029 

9.196 

8.772 

Each  slope  may  be  obtained  by  multiplying  the  corresponding 
value  in  the  table  by  0.003225- 
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TAB1£  V 


MEASURED  SLOPES  IN  PITCH 


sr 


OF  NA'nJRAL-VIBRATION-MQDE  SHAPES® 


i 

S 

Semi  span 
station, 
percent 

r  =  1 

(25  percent 
chord) 

r  »  2 

(50  percent 
chord) 

r  =  3 

(75  percent 
chord) 

r  a  4 
( trailing 
edge) 

1 

20 

0.1063 

0.204 

0.0968 

-0.166 

2 

35 

.572 

.540 

.484 

.0445 

5 

50 

1.19 

1.16 

.854 

.610 

1 

4 

65 

2.015 

1.83 

1.655 

l.4l6 

5 

80 

2.735 

2.69 

2.62 

2.47 

6 

90 

3.61 

3.355 

3.54 

3.315 

1 

20 

-0.205 

-0.221 

-0.174 

-0.164 

2 

35 

-.446 

-.2895 

-.^73 

-.755 

3 

50 

.1208 

-.0836 

-.348 

-1.13 

2 

4 

65 

1.02 

.931 

.^99 

-.2455 

5 

80 

3.5J+ 

3.01 

2.98 

2.32 

6 

90 

6.18 

5.80 

5.59 

5.^9 

1 

20 

-0.074 

-1.29 

-1.295 

-1.78 

2 

35 

-.556 

-1.713 

-3.05 

-4.14 

3 

50 

-2.015 

-3.^95 

-5.05 

-6.50 

J 

4 

65 

-4.19 

-5.26 

-6.65 

-7.69 

5 

80 

-6.26 

-6.94 

-7.70 

-8.15 

6 

90 

-7.55 

-7.58 

-7.85 

-7.98 

1 

20 

-0.172 

0.067 

0.335 

0.835 

2 

35 

.028 

.084 

.483 

1.426 

5 

50 

-.474 

-.301 

.082 

.479 

4 

4 

65 

.791 

-.483 

-.7^5 

-1.270 

5 

80 

.790 

.767 

.371 

-.461 

6 

^  1 

4.515 

4.172 

3.640 

2.937 

a 

Each  slope  may  be  obtained  by  multiplying  the  corresponding 
value  in  the  table  by  0.003225 • 


DATA  04PLOYED  IN  NUMERICAL  INTEGRATION 
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Control  points  are  shown  in  figure  l(a). 
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chord  chord 


(a)  Control  points. 


(b)  Numerical  integration. 

Figure  1.-  Lattices  used  in  evaluation  of  elements  of  flutter  detezmi* 

nant  for  delta  semispan  ving. 


23.35 
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Figure  2.-  Construction  details  of  delta-wing  model 


HhCA  tn  1+595 
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•w. 

Ib/sq  in. 


(a)  Chordvlse  distribution. 


Figure  5*- 


Weight  distributions  of  model 


z.* 
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Figure  4.-  Continued 
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Figure  7.-  Schematic  diagram  of  appeuratus  used  to  measure  natural-vibrat ion-mode  shapes.  All 
directions  are  considered  positive  as  shown;  z  is  considered  positive  downward. 


